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Most investigations of multipartite entanglement have been concerned with 
temporal modes of the electromagnetic field, and have neglected its spatial 
structure. We present a simple model which allows to generate tripartite en- 
tanglement between spatial modes by parametric down-conversion with two 
symmetrically-tilted plane waves serving as a pump. The characteristics of 
this entanglement are investigated. We also discuss the generalization of our 
scheme to 2N + 1-partite entanglement using 2N symmetrically-tilted plane 
pump waves. Another interesting feature is the possibility of entanglement 
localization in just two spatial modes. 
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1. Introduction 

Multipartite entanglement, as the name suggests, is the entanglement between more than two parties. The existence 
of this type of entanglement for both discrete and continuous quantum variables is one of the most striking predictions 
of quantum mechanics and is very unusual from the point of view of classical physics. It is therefore not surprising 
that nowadays multipartite entanglement is a very active research area of quantum optics and quantum information. 
The interest to multipartite entanglement is motivated not only by its fundamental nature but also by its potential 
for applications in quantum communication technologies. 

Numerous theoretical publications have addressed the characterization of multipartite entanglement. Coffman, 
Kundu and Wootters [1] have established for a three-qubit system and conjectured for A^-qubit systems the so-called 
monogamy of quantum entanglement, constraining the maximum entanglement between partitions of a multiparty 
system. Later on, Adesso and lUuminati [2]- [5] have introduced the continuous- variable tangle as a measure of 
multipartite entanglement for continuous variable (CV) multimode Gaussian states and have demonstrated that it 
satisfies the Coffman-Kundu- Wootters monogamy inequality. The conjecture of Ref. [1] has been proven by Osborne 
and Verstraete [6] . The corresponding proof for Gaussian states is in Hiroshima, Adesso and Illuminati [7] . 

Among the potential applications of CV multipartite entanglement let us mention quantum teleportation networks, 
quantum telecloning, controlled quantum dense coding, and quantum secret sharing [8]. 

Several generation schemes for CV multipartite entanglement have been proposed in the literature theoretically 
and realized experimentally in recent years. Among the first ones was the "passive" optical scheme using squeezed 
states mixed with beam splitters [9]. Then came the "active" schemes in which multipartite entanglement is created 
as a result of parametric interaction of several optical waves such as cascaded/concurrent [10]-[12], interlinked [13], 
[14] or consecutive parametric interactions [15]. All of these schemes for generation of CV multipartite entanglement 
are considering the temporal modes of the electromagnetic field and neglect its spatial structure. The question arises 
on whether the spatial modes can also serve for the creation of CV multipartite entanglement [16]- [19]. In Ref. [3], a 
physical implementation is given in terms of 2A^-1 beam splitters and 2 N single- mode squeezed input states, itself 
based on a previous work [20]. 

In this paper we propose a simple scheme which corresponds to the "active" creation of tripartite entanglement be- 
tween spatial modes of the electromagnetic field in the process of parametric interaction. We create spatial tripartite 
entanglement by pumping a nonlinear parametric medium by a coherent combination of several tilted plane monochro- 
matic waves, what we call a spatially-structured pump. Since the pump photon can be extracted from any of these 
waves and the pair of down-converted photons emitted in different directions according to the phase-matching condi- 
tion, our scheme allows for the creation of tripartite entanglement between spatial modes of the down-converted field. 
An interesting feature of our scheme is the possibility of localizing the created spatial tripartite entanglement in just 
two well-defined spatial modes formed as a linear combination of all the modes participating in the down-conversion 
process. 

In Sec. II, we shall consider explicitly the case of tripartite entanglement which is produced by parametric down- 
conversion with two symmetrically-tilted plane waves serving as a pump. The genuine tripartite entanglement is 
analyzed in subsection A and its localization on two modes is studied in subsection B. We then discuss in Sec. Ill 
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Fig. 1. (Color online) Scheme for generating tripartite entanglement between spatial modes. The tilted pumps have 
wave vectors kp(±go)- The transverse (vertical) components are ±qo. 

a generalization of our scheme to 2N + 1-partite entanglement using 2N symmetrically-tilted pump waves and draw 
some conclusions. 

2. Spatial tripartite entanglement 

A. Parametric down- conversion with two symmetrically-tilted plane waves 
The scheme of parametric down-conversion used here is shown in Fig. 1. 

We consider two plane pump waves with wave vectors kp(±<7o) where ±(70 ai'^ the pertaining transverse components 
(which are vertical on Fig. 1) at the input to the crystal. The pump field reads then 



'p(a;) =a(e**"-f e-**"), (1) 



with a the amplitude of the pump, chosen as a real number in our case. We shall denote the corresponding spatial 
Fourier components of the corresponding field operators as follows 

a{z,q) EE aoiz) 

d{z,±qo + q) = a±{z). (2) 

These modes are depicted in Fig. 1: the long (blue) arrows pertains to a±{0), the short one to aj(0). 

Due to the parametric interaction in the nonlinear medium these operators are coupled and obey the following 
equations in the rotating wave approximation [21]: 



d 

— ao = a 
az 

d_ 

dz 



\(d\ + a^_Y'^'- , 



—a+ = aXaUe 



= aXale'^', (3) 
dz 

where A is the coupling constant and A is the phase mismatch. 

The solution to this set of equations for the fields at the output of the crystal, z = I, can be written analytically as 
follows 

= i[C/(0 + l]a+(0) + i[C/(0-l]a-(0) + ^n04(0), 
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U{l)an{Q) + ^V{l)[a}+(Q) 



l]a-(0) + -^F(04(0), 



(4) 



with the functions 



U{z) 
V{z) 



JAz/2 



cosh(rz) 



*^sinh(r2;; 



e'^^/^^sinh(rz), 



(5) 



where F 



We shall be concerned for simplicity with the case of zero phase mismatch, A = 0, and 
will be considering only the zeroth spatial Fourier components of the field, q = Q. Physically, this corresponds to 
photodetection of the light field by a single large photodetector without spatial resolution. In this case, we have 
U{1) = coshr and V{1) = sinhr with r = y/2aXL being the squeezing parameter. 

When the input states of the three interacting modes are vacuum states, the output states remain Gaussian. Hence, 
they are completely described by the output covariance matrix a associated with the following quadrature components 

xo(^) = 2Reao(^) 

Po{t) = 2Imao(^) 

x±{l) EE 2Rea±(^) 

p±{l) = 2Ima±(^). (6) 



Its matrix elements read 



Tr 



(Ae.AC, + A^-Ae.) /2 



(7) 



where p is the output density matrix and A^.^ = £,i — {^dj with the i-th component of the vector ^ — 

{xo{£),po{£), x-^-{£),p-^-{£), x^{£),p^{£)). Taking the solution Q into account, the covariance matrix at the output of 
the crystal reads explicitly 



f cosh(2r) 
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(8) 



Notice that this matrix is bisymmetric as it is invariant under the permutation of the quadratures pertaining to the 
modes + and -: x+{£),p+{£) O X-{£),p-{£) . Hence, it has the tripartite entanglement structure of covariance 
matrices associated to bisymmetric (l+2)-mode Gaussian states considered in Ref. [3]. Actually, we can show that 
the output state obtained here exhibits genuine tripartite entanglement in the sense of Ref. [20]. For that purpose, 
we have to verify that the following condition on covariance matrix elements is violated 



C 



xoi£) 




Mi) 



P-{£)+P-{£) 
V2 




From dS]) one deduces that 

C = 4{cosh(2r) - sinh(2r)} = 4e 
which is smaller than 1/2 if the squeezing parameter r > 



-2r 



(9) 



(10) 

y/2a\l therefore entails that 



In 2. Recalling that r 

the output state p produced by the above parametric down-conversion process with two symmetrically-tilted plane 
waves exhibits genuine tripartite entanglement when the pump amplitude a, the coupling parameter A and the crystal 
length £ are such that 

a\£ > w 0.735. (11) 
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B. Two-mode entanglement localization 



We shall show that the tripartite entanglement generated by the above process can be concentrated on two modes by 
local unitary transformations as was discovered by Serafini, Adesso and Illuminati [3]. Let us consider the following 
quadratures 



4W 



(12) 



V2 



72 



which are obtained from the previous ones by a unitary transformation corresponding to the action of a beamsplitter. 
Accordingly, the output covariance matrix a is thus transformed by congruence, 



S^aS, 



(13) 



where S is the following symplectic transformation. 
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(14) 
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with 1\ the 2x2 identity matrix. The new covariance matrix reads then 

/ cosh(2r) sinh(2r) 

, _ I cosh(2r) 

^ ^ sinh(2r) cosh(2r) 

V -sinh(2r) 




(15) 



The separability of the pertaining state p' can be determined by considering a partition of the system in two 
subsystems and investigating the positivity of the partially transposed matrix p\ obtained upon transposing the 
variables of only one of the two subsystems. The positivity of partial transposition (PPT) is a necessary condition 
for the separability of any bipartite quantum state [22], [23]. It is also sufficient for the separability of (1 + n)— mode 
Gaussian states [24], [25]. 

The covariance matrix a' of the partially transposed state p' with respect to one subsystem is obtained [24] by 
changing the signs of the quadratures belonging to that subsystem. Here, because of the block-diagonal structure 
of a' we can restrict the analysis to the first block and consider the transposition of the mode 1, i. e. change the sign 
of pi, 



cosh(2r) sinh(2r) 

cosh(2r) sinh(2r) 

sinh(2r) cosh(2r) 

sinh(2r) cosh(2r) 



(16) 



The symplectic eigenvalues {vi\ of the covariance matrix a' are v± — e and 1. The necessary and sufficient PPT 
condition for the separability of the state p' amounts to having > 1 for all i. We can therefore focus on the smallest 
eigenvalue v-. The extent to which the separability of the state is violated, is measured by £js/{p), the logarithmic 
negativity of p' , defined as the logarithm of the trace norm of p\ 

£aa(p') = ln||/5'||i =max(0,-lnz>_). (17) 

Using the calculated symplectic eigenvalues {vi} we obtain a very simple expression for the logarithmic negativity, 

£m{p') - 2r. (18) 



This implies that the central mode is entangled with the mode 1 which is precisely the combination of modes + and 
- that was used in the condition ^ in order to show the existence of genuine tripartite entanglement. 
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Fig. 2. (Color online) Generalization of the scheme for 2N symmetrically tilted pump waves with N = 2 and 4. The 
little circles represent the projections of the pump wave vectors in the xy plane of the crystal entering face. 



3. Discussion and conclusion 

Our scheme can be generalized to multipartite entanglement localization using 2N symmetrically tilted plane waves 
as pump modes in the parametric down-conversion process. An illustration is provided in Fig. 2 for the case N — 2 
and 4 which correspond to 4 and 8 pump waves. 



The projections of the pump wave vectors in the xy plane of the crystal entering face are depicted as little circles. In 
this case, the system contains 2N + 1 coupled modes which are described by equations similar to 

We can solve them in the rotating wave approximation and investigate the properties of the multipartite entangle- 
ment in this case (the details will be provided in a forthcoming publication) . The covariance matrix a at the output of 
the crystal is again bisymmetric as those considered in Ref . [3] . By an orthogonal transformation of the field operators 
which generalizes the one achieved in the tripartite case, the covariance matrix cr can be brought to a form similar to 
(|15p . namely the tensor product of 

a 4 X 4 covariance matrix and a 2{2N — 1) x 2{2N — 1) identity matrix. The 4x4 covariance matrix has exactly 
the same form as the one in ([8]) with r replaced by ^/Nt. Hence, it corresponds to that of two entangled modes with 
the squeezing parameter enhanced by n/TV and 2N — 1 modes in the vacuum state. For the logarithmic negativity we 
indeed obtain £_\f{p') = 2\/Nr. Physically, this means that by mixing different spatial modes with beam splitters, 
we can "localize" the entanglement distributed initially among all the 27V -f 1 spatial modes in only two well-defined 
modes, formed by the linear combinations of the initial ones. Interestingly, this entanglement localization results in 
an enhancement of the entanglement by a factor V^- 

To conclude, our "active" optical scheme relies on the use of a spatially-structured pumping of the parametric 
medium. It provides a realistic proposal for the experimental realization of multipartite entanglement as well as its 
localization. The possibility of entanglement localization was predicted mathematically in Ref. [3] and a physical 
implementation was also given in terms of 2N-1 beam splitters and 2N single- mode squeezed inputs [20]. We have 
also established the scaling behavior of the logarithmic negativity with the number of modes. Note that due to the 
asymmetric coupling between the modes in our scheme, this scaling is a relatively slowly growing function of the 
number of modes. An interesting question is to determine the ultimate rate at which the logarithmic negativity may 
grow with the number of modes. This subject will be addressed elsewhere. 
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